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^ I For a semialgebraic set K in E n , let Pd(K) = {/ G R[a;]<d : /(it) > OVit £ if} be the 

' cone of polynomials in x G M™ of degrees at most e? that are nonnegative on K . This paper 

studies the geometry of its boundary dPd(K). When K = W 1 and d is even, we show 
that its boundary dPd(K) lies on the irreducible hypcrsurface defined by the discriminant 
A(/) of /. When K = {x G R" : gi(x) = ■ ■ ■ = g m {x) = 0} is a real algebraic variety, wc 
U ' show that dPd(K) lies on the hypersurface defined by the discriminant A(/, gi, . . . ,g m ) 

^ <—| of /, gi, . . . ,g m . When K is a general semialgebraic set, we show that dPd(K) lies on a 

union of hypersurfaces defined by the discriminantal equations. Explicit formulae for the 
degrees of these hypersurfaces and discriminants are given. We also prove that typically 
Pd{K) does not have a barrier of type — log ip(f) when ip(f) is required to be a polynomial, 
but such a barrier exits if tp(f) is allowed to be semialgebraic. Some illustrating examples 
04 ' are shown. 
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1 Introduction 

Let K be a semialgebraic set in R n , and Pd(K) be the cone of multivariate polynomials in 
x 6 M. n that are nonnegative on K and have degrees at most d, that is, 



Pd{K) = {/ G R[x]< d : f(u) > V it G K} . 



A very natural question is what is the boundary of Pd(K)? What kind of equation does it 
satisfy? Can we find a nice barrier function for Pd{K)l This paper discusses these issues. 

A polynomial f(x) in x G M. n is said to be nonnegative or positive semidefinite (psd) on K 
if the evaluation f(u) > for every u G K. When K = W 1 and d is even, an f(x) G Pd(W l ) 
is called a nonnegative polynomial or psd polynomial. When K = W\_ is the nonnegative 
orthant, an f(x) G P^MIM is called a co-positive polynomial. Typically, it is quite difficult 
to check the membership of the cone Pd{K). In case of K = M. n , for any even d > 2, it is 
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NP-hard to check the membership of Pd(W"). In case of K = W], for any d > 1, it is NP-hard 
to check the membership of In practical applications, people usually do not check 

the membership of Pd{K) directly, and instead check sufficient conditions like sum of square 
(SOS) type representations (a polynomial is SOS if it is a finite summation of squares of 
other polynomials). There is much work on applying SOS type certificates to approximate 
the cone P d {K). We refer to p2J HI [20j EH E2J [27] . However, there is relatively few work 
on studying the cone Pd(K) and its boundary dPd(K) directly. The geometric properties of 
dPd(K) are known very little. 

When K = W 1 and d = 2, P2(W l ) reduces to the cone of positive semidefinite matrices, 
because a quadratic polynomial f(x) is nonnegative everywhere if and only if its associated 
symmetric matrix A y 0(positive semidefinite). The boundary of P2(M. n ) consists of / whose 
corresponding A is positive semidefinite and singular, which lies on the irreducible determi- 
nantal hypersurface det(A) = 0. Its degree is equal to the length of matrix A. A typical 
barrier function for P2(W l ) is — logdet(j4). Note that det(A) is a polynomial in the coeffi- 
cients of f(x). Do we have a similar result for Pd(K) when K ^ W 1 or d > 2? Clearly, 
when K = W 1 and d > 2, we need to generalize the definition of determinants for quadratic 
polynomials to higher degree polynomials. There has been classical work in this area like 
[6]. The "determinants" for polynomials of degree 3 or higher are called discriminants. The 
discriminant A(/) of a single homogeneous polynomial (also called form) f(x) is defined such 
that A(/) = if and only if f(x) has a nonzero critical point. For a general semialgebraic set 
K, to study dPd(K), we need to define the discriminant A(/o, . . . , f m ) of several polynomials 
/0) • • • j fm- As we will see in this paper, the discriminant plays a fundamental role in studying 
P d (K). 

Recently, there are arising interests in the new area of convex algebraic geometry. The 
geometry of convex (also including nonconvex) optimization problems would be studied by 
using algebraic methods. There is much work in this field, like maximum likelihood esti- 
mation [2], /c-ellipse [IT], semidefinite programming \19\ 123]. matrix cubes [18] . polynomial 
optimization |16j . statistical models and matrix completion [29], convex hulls [101 [24] 126] . 
theta bodies [TJ. In this paper, we study the geometry of the cone Pd{K) by using algebraic 
methods, and find its new properties. 

Contributions The cone Pd(K) is a semialgebraic set, and its boundary dPd(K) is a 
hypersurface defined by a certain polynomial equation. To study this hypersurface, we 
need to define the discriminant A(/o, ■ ■ ■ , fm) for several forms /o, . . . , / m , which satisfies 
A(/o, . . . , f m ) = if and only if fo(x) = • • • = f m (x) = has a nonzero singular solution. 
This will be shown in Section [3] When K = W 1 and d > 2 is even, we prove that dPd(M. n ) lies 
on the irreducible discriminantal hypersurface A(/) = 0, which will be shown in Section [5J 
When K = {x G M n : g±(x) = ■ ■ ■ = g m (x) = 0} is a real algebraic variety, we show that 
dPd(K) lies on the discriminantal hypersurface A(/oj ■ ■ ■ , fm) = 0, which will be shown in 
Section [5] When K is a general semialgebraic set, we show that dPd(K) lies on a union of 
several discriminantal hyper surf aces, which will be shown in Section [6] Explicit formulae for 
the degrees of these hypersurfaces will also be shown. Generally, we show that Pd(K) does 
not have a barrier of type — log <p(f) when (p(f) is required to be a polynomial, but such a 
barrier exits if f(f) is allowed to be semialgebraic. For the convenience of readers, we include 
some preliminaries about elementary algebraic geometry, discriminants and resultants. This 
will be shown in Section [2] 
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2 Some preliminaries 



2.1 Notations 

The symbol N (resp., M) denotes the set of nonnegative integers (resp., real numbers), and 
denotes the nonnegative orthant of M n . For integer n > 0, [n] denotes the set {1, . . . , n}. 
For x G R n , Xj denotes the i-th component of x, that is, x = (xi, . . . ,x n ), and x denotes 

(xq, xx,..., x n ). For a G N n , denote |a| = ax H h an- For i£l n and a G N n , x a denotes 

x" 1 • • -x" n . The [x rf ] denotes the column vector of all monomials of degree d, i.e., [x d ] T = 

[xf x(~ 1 x 2 x*]. The symbol R[x] = R[xx,...,x n ] (resp. C[x] = C[xx, . . . ,x n ]) 

denotes the ring of polynomials in {x\, . . . ,x n ) with real (resp. complex) coefficients; R[x] = 
R[xo,xi, . . . ,x n ] and C[x] = C[xq,xi, . . . ,x n ] are defined similarly. A polynomial is called 
a form if it is homogeneous. The K[x]d (resp. M[x]d) denotes the subspace of homogeneous 
polynomials in R[x] (resp. R[x]) of degree d, and M[x]<^ = M[x]o + R[x]x + • • • + R[x]d- For 
a polynomial f(x) of degree d, f h {x) denotes its homogenization Xq/(x/xo). For a tuple 
g = (gx, ■ ■ ■ ,g m ) of polynomials, denote g h = (g^, . . . For a finite set S, \S\ denotes 

its cardinality. For a general set S C R n , int(S) denotes its interior, and dS denotes its 
boundary in standard Euclidean topology. For a matrix A, A T denotes its transpose. For a 
symmetric matrix X, X y (resp., X y 0) means X is positive semidefinite (resp. positive 
definite). For u G M. N , \\u\\2 = V u T u denotes the standard Euclidean norm. 

2.2 Ideals and varieties 

In this subsection we give a brief review about ideals and varieties in elementary algebraic 
geometry. We refer to [3j [9] for more details. 

A subset / of C[x] is called an ideal if p ■ q G / for any p G M[x] and q £ I. For 
gx, ■ ■ ■ ,g m G C[x], (gx, • • • , g m ) denotes the smallest ideal containing every gi. The gx,...,g m 
are called generators of (g\ , • • • , g m ), or equivalently, (gx, ■ ■ ■ , g m ) is generated by gx, ■ ■ ■ ,g m - 
Every ideal in C[x] is generated by a finite number of polynomials. 

An algebraic variety is a subset of C n that is defined by a finite set of polynomial equations. 
Sometimes, an algebraic variety is just called a variety. Let g = (gx, ■ ■ ■ ,g m ) be a tuple of 
polynomials in R[x]. Define 

V{jg) = {x G C n : g\(x) = ■ ■ ■ = g m (x) = 0}. 

In optimization, we are more interested in real solutions. Define 

Vn(g) = {x G R n : gx(x) = ■ ■ ■ = g m (x) = 0}. 

It is called a real algebraic variety. Clearly, Vu(g) C V(g). If I = (gx, ■ ■ ■ ,g m ), we define 
V(I) = V(g). 

Given V C C n , the set of all polynomials vanishing on V is an ideal and denoted by 

I(V) = {h G C[x] : h(u) = V u G V}. 
Clearly, if V = V(I) and p G /, then p G i"(V). The following is a reverse to this fact. 

Theorem 2.1 (Hilbert's Nullstellensatz). Let I C C[x] be an ideal. If p G I(V), then p k G J 
for some integer k > 0. 
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Given a subset S C C n , the smallest variety V C C n containing S is called the Zariski 
closure of S, and is denoted by Zar(S). For instance, for S = {x G M 2 : x\ + x\ = 1, xi > 
0,^2 > 0}, its Zariski closure is the variety {x G C 2 : x\ + x\ = 1}. In the Zariski topology 
on C n , the varieties are closed sets, and the complements of varieties are open sets. 

The varieties in the above are also called affine varieties, because they are defined in the 
vector space C n or W 1 . We also need projective varieties that are often more convenient in 
algebraic geometry. Let P n be the n-dimensional complex projective space, where each point 
x G P n is a family of nonzero vectors x = {xq,x\, . . . ,x n ) that are parallel to each other. 
A set U in P" is called a projective variety if it is defined by finitely many homogeneous 
polynomial equations. For given forms pi(x), . . . ,p m (x), denote the projective variety 

Vwipu ■ ■ ■ ,Pm) = {x G P n : pi(x) = ■■■= p r (x) = 0} . 

In particular, if m = 1, Vp(p\) is called a hypersurface. Furthermore, if p\ has degree one, 
Vp(pi) is called a hyperplane. In the Zariski topology on P n , the projective varieties are 
closed sets, and their complements are open sets. 

A variety V is irreducible if there exist no proper subvarieties V\,V2 of V such that 
V = V\ U V2- The dimension of a variety U is the biggest integer t such that U = Uq D U\ D 
• • • D Ug where every Ui is an irreducible variety. For an ideal / C C[x], its dimension is 
defined to be the dimension of its variety V(I). It is zero-dimensional if and only if V(I) is 
finite. 

Let V be a projective variety of dimension I in P n and /(V) = . . . , f r ). The singular 
locus V S i n g is defined to be the variety 

V sing = {w £V : rank J(/ x , . . . , f r ) < n - £ at w} , 

where J(/i, . . . , f r ) denotes the Jacobian matrix of fx, . . . , f r . The points in V s i ng are called 
singular points of V. If V s % n g = 0, we say V is smooth. When V is an affine variety, its 
singular locus and singular points are defined similarly. 

2.3 Discriminants and resultants 

In this subsection, we review some basics about discriminants and resultants for multivariate 
polynomials. We refer to [6] for more details. 

Let f(x) be a polynomial in x = (x\, . . . ,x n ) and u G C n be a complex zero point of 
f(x), i.e., f(u) = 0. We say u is a critical zero of / if V x f(u) = 0. Not every polynomial 
has a critical complex zero. In the univariate case, if f(x) = ax 2 + bx + c is quadratic 
and has a critical complex zero, then its discriminant b 2 — 4ac = 0. In the multivariate 
case, if f(x) = x T Ax is quadratic and A is symmetric, then f(x) has a nonzero complex 
critical point if and only if its determinant det(^4) = 0. The above can be generalized 
to polynomials of higher degrees. In [6], the discriminants have been defined for general 
multivariate polynomials. 

For convenience, let f(x) be a form in x = (x±, . . . , x n ). The discriminant A(/) is a 
polynomial in the coefficients of / satisfying 

A(/) =0 ^ 3 u G C n \{0} : V/(u) = 0. 
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The discriminant A(/) is homogeneous, irreducible and has integer coefficients. It is unique 
up to a sign if all its integer coefficients are coprime. When deg(/) = d, A(/) has degree 
n(d — l) n_1 . For instance, when n = 2 and d = 3, we have the formula (see [6J Chap. 12]) 



A(axf + bx\x 2 + cx x x\ + dx\) = b 2 c 2 - 4ac 3 - 46 3 d + 18abcd - 21a 2 d 2 . 

A more general definition than discriminant is resultant. Let f\, . . . , f n be forms in x S M. n . 
The resultant Res(f\, . . . , f n ) is a polynomial in the coefficients of fi, ■ ■ ■ , f n satisfying 

Res(f t,..., f n ) = ^ 3u G C n \{0} : /i(u) = ••• = /„ (u) = 0. 

The polynomial Res(f\, . . . , /„) is homogeneous, irreducible and has integer coefficients. It is 
unique up to a sign if all its coefficients are coprime. If fi has degree di, then Res(f±, . . . , / n ) 
is homogeneous in every of degree d\ ■ ■ ■ • • • d n , and its total degree is 

d x ■ ■ ■ d n (d^ 1 + ■■■ + d- 1 ) . 

In case of n = 2, a general formula for Res(fi, . . . , f n ) is given in [281 Sec. 4.1]. For instance, 
if fi(x) = ax\ + bx\X2 + cx\ and f2(x) = dx\ + exiX2 + fx 2 , then 

-Res(/i, / 2 ) = c 2 d 2 - bcde + ace 2 + 6 2 oy - 2acd/ - abef + a 2 f 2 . 

We would like to to remark that the discriminant is a specialization of resultant. A form 
f(x) has a nonzero complex critical point if and only if 

df(x) = = dfjx) = Q 
dx\ dx n 

has a nonzero complex solution. So A(/) = r] ■ Res(-7^, . . . , for a scalar -q ^ 0. 

In many situations, we often handle nonhomogeneous polynomials. The discriminants 
and resultants would also be defined for them. Let f{x) be a general polynomial in x = 
(x\, . . . , x n ), and the form f h (x) in x = (xo, X\, . . . , x n ) be its homogenization. The dis- 
criminant A(/) of f(x) is then defined to be A(f h ) . Observe that if u £ C" is a critical 
zero point of /, i.e., f(u) = and V x f(u) = 0, then we must have Vxf h {u) = 0. Here 
u = (1, • • • ,u n ). To see this point, recall the Euler's formula (suppose deg(/) = d) 

OXq oxi ox n 

Since f h (u) = f(u),U^ = ^ , . . . , *0 = it holds that V,f\u) = 0. It is 

possible that A(/) = while / does not have a critical zero point, because Vxf h (x) = 
might have a solution at infinity xq = 0. 

The resultants are similarly defined for nonhomogeneous polynomials. Let /o, fx, ■ ■ ■ , f n 
be general polynomials in x = (x±, . . . , x n ). The resultant i?es(/o, f±, . . . , f n ) is then defined 
to be Res{fQ, f±, . . . , f%). Here each form f^{x) is the homogenization of fi{x). Clearly, if 
the polynomial system 

/o(x) = /i(x) = --- = /«(a;) = 
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has a solution in C n , then the homogeneous system 

f^x) = f 1 h (x) = --- = f^x) = 

has a solution in P n . The reverse is not always true, because the latter might have a solution 
at infinity xq = 0. 

There are systemic procedures to compute resultants (hence including discriminants) for 
general polynomials. We refer to [H Chap. 3], [6j Sec. 4, Chap. 3], and [28j Chap. 4]. 



3 Discriminants of several polynomials 

In this section, we assume fo(x), fi(x), . . . , f m (x) are forms in x = (xq, x\, . . . , x n ) of degrees 
do, d\, . . . , d m respectively, and m < n. Denote / = (fo, /i, . . . , / m ). If every /j has generic 
coefficients, the polynomial system 

fo(x) = ■ ■ ■ = f m (S) = (3.1) 

has no singular solution in P n , that is, for any ti£P n satisfying (|3.ip . the Jacobian 

J,(u) := [V 4 / (iQ Vx/i(fi) ••• Vi/ m (u)] 

has full rank. For some particular /, (|3.ip might have a singular solution. Define 



W(d ,...,d r , 



3u 6 P n s.i 

(/0, • • • Jm) € n R [ f ]^ : /o(«) = • • • = /»»»(«) = 
i=o rankJj(u) < m 



When every di = 1, W(l, . . . , 1) consists of all vector tuples (/o, . . . , f m ) such that fo, ■ ■ ■ , f m 
are linearly dependent. Thus W(l, ■■-,!) consists of all (n + 1) x (m + 1) matrices whose 
ranks are at most m, which is a determinantal variety of codimension n + 1 — m. It is not 
a hypersurface when m < n — 1. When every = d > 1, W(d, . . . ,d) consists of all tuples 
(fo, ■ ■ ■ i fm) such that the multi-homogeneous form in (i, A) (here A = (Ao, Ai, . . . , A m )) 

C(x, A) := X fo(x) + \ifi(x) H h X m f m (x) 

has a critical point in the product of projective spaces P n x P m . As is known, the multi- 
homogeneous form C(x, A) has a critical point in P n x P m if and only if its discriminant 
vanishes (see (HJ Section 2B, Chap. 13]). So W(d, . . . , d) is a hypersurface. When the d,'s are 
not equal and at least one di > 1, W(do, . . . , d m ) is also a hypersurface, which is a consequence 
of Theorem 4.8 of Looijenga [14]. This fact was kindly pointed out to the author by Kristian 
Ranestad. So we assume at least one di > 1, and then W(do, . . . ,d rn ) is a hypersurface. 
Let A(/o, fi, . . . , f m ) be a defining polynomial of the lowest degree for W(do, . . . , d m ). It is 
unique up to a constant factor and satisfies 

(f ,...,f m )eW(d ,...,d m ) A(/ , /!,..., / m ) = 0. (3.2) 

For convenience, we also call A(/o, fx, . . . , f m ) the discriminant of forms fo(x), . . . , f m (x). 
When m = 0, A(/o, f\, . . . , fm) becomes the standard discriminant of a single form, which 



G 



has degree (n + l)(do — l) n . So A(/o, /i, . . . , f m ) can be thought of as a generalization 
of A(/o). In the rest of this section, we are going to prove a general degree formula for 
A(/o> /i, ■ ■ ■ , / m ). 

For every integer k > 0, denote by the fe-th complete symmetric polynomial 

Sfc(ai, . . . = ^ a 4 1 1 ---aj t . 

iiH Mt=fc 

Let iif(x) G M[x]( ri+1 W m+1 ) be a matrix polynomial such that its every entry Hij(x) is 
homogeneous and all the entries of its every column have the same degree. Define 

V m {H) = {x e P n : ranktf(x) < m}. (3.3) 

Theorem 3.1. Suppose every di > 0, at least one d{ > 1, and m < n. Then the discriminant 
A(/o, . . . , f m ) has the following properties: 

a) For every k = 0, . . . , m, A(fo, fi, ■ ■ ■ , fm) is homogeneous in fk- It also holds that 

A(/o, ■ ■ • , fm) = whenever fi = f j for i / j. 

b) For every k = 0, . . . , m, the degree of A(/ , /i, . . . , f m ) in fk is 

do ■ ■ ■ d k ■ ■ ■ d m ■ S n -m (do - 1, . . . , d k - 1, . . . , d m - lj ■ (3.4) 

In the above, dk means dk is missing, and a means a is repeated twice. Thus the total 
degree of A(/ , /i, . . . , f m ) is 




c) For fixed fi, . . . , f m , A(/o, /i, . . . , f m ) is identically zero in fo if and only if the projec- 
tive variety Vp(/i, . . . , f m ) has a positive dimensional singular locus. 

Proof, a) Note that for any scalar a ^ 0, (fo, • • • , f m ) £ W(do, . . . , d m ) if and only if 

(fo, ■ - - , fk-l,Ctfk, fk+l, ■ ■ fm) G W(d , ■ ■ ■ ,d m )- 

So, by relation (13. 2ft . A(/o, . . . , f m ) must be homogeneous in every fk- 

If /; = fj for some distinct i,j, say i = 0,j = 1, then (fo,...,f m ) G VF(do, . . . , d m ) 
because the polynomial system (|3.ip must have a solution in P n (it has only m — 1 < n 
distinct equations) and its Jacobian is singular (its first two columns are same). 

b) For convenience, we only prove the degree formula for k = 0. Choose generic forms 
fo, ■ ■ ■ , fm of degrees do, ■ ■ ■ ,d m respectively, and another generic form h of degree do- Then 
the degree of A(/o, f\, . . . , f m ) in fo is equal to the number of scalars 7 such that 

A(/ + 7 h, /!,..., / m ) = 0. (3.6) 

Since the /,'s are generic, A(/i, . . . , f m ) / and hence Mp(/i 5 • • • , fm) is nonsingular. 
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Lemma 3.2. The condition \3. 6\) is equivalent to 

3n£P n ,3 7 GC: h{u) = ■ ■ ■ = f m (u) = 0, 
rank [V x fo(u) + 7 V ' x h{u) ^xfi(u) ■ ■ • V £ / m (u)] < m. 

Furthermore, every u satisfying \3. 7\ ) determines 7 uniquely. 

Proof. By relation ()3.2p . (13. 6p clearly implies (13. 7p . So we only prove the reverse. Sup- 
pose (|3.7p is satisfied by some u and 7. The rank condition in (|3.7p implies there exists 
(/i , Ml, • • • , Mm) 7^ satisfying 

Mo(Vx/ (ii) + 7Vi/i(u)) + H\Vxh{u) H h /i m V i / m (u) = 0. 

Since Vp(/i, . . . , / m ) is nonsingular, we must have fj,Q / and can scale /^o = 1- By Euler's 
formula (|2.ip . premultiplying u T in the above equation gives 

do(fo(u) + 7/1(1*)) + Hidih(u) H h ^ m d m f m (u) = 0. 

Thus the equations in (|3.7p imply /(u) + 7/1(71) = 0, and thus (I3.6j) holds by relation (13. 2p . 

Now we prove each u in (|3.7p uniquely determines 7. If /i(u) / 0, we know 7 = —f(u)/h(u) 
from the above. If h{u) = 0, because Vp(/i, /1, . . . , / m ) is nonsingular (/i and fi are all generic), 
we can generally assume the first m + 1 rows of the Jacobian of h, fi, ■ ■ ■ , f m at u are linearly 
independent, which is denoted by [b F] with b G (C m + 1 and F G O m+1 ' xm . Denote by a 
the first m + 1 entries of V^/o(w). Then, det [6 F] 7^ and ()3.7p implies 

det [a + 76 F] = det [a F] + 7 det [6 F] = 0. 

So 7 = — det [a F] / det [b F] . There is a unique 7 for every u in (|3.7p . □ 
Clearly, (13.7j) is also equivalent to 

BuGP": / 1 (u) = ... = / ro ( u ) = 0, 
rank [Vxfo(u) V x h(u) S/ X fi(u) ••• V x f m (u)] < m + 1. 

Let J be the Jacobian matrix in the above. By Lemma 13.21 the degree of A(/o, /1, . . . , f m ) 
in /o is equal to the cardinality of 

u-.= v m+1 (j)nv w (f f m ). 

The variety Vp(/i, . . . , f m ) is smooth, has codimension m and degree d\ ■ ■ ■ d m . Since every 
fi and h are generic, XX m+ i(J) is also smooth, has dimension m and intersects Vp(/i, . . . , / m ) 
transversely. So C/ is a finite variety. We refer to Proposition 2.1 and Theorem 2.2 in |16j 
for more details about this fact. The degree of the determinantal variety P m +i(J) is (cf. 
Proposition A. 6 of [16]) 

S n -m{do — 1, ^o — 1) d\ — 1, . . . , d m — 1). 

By Bezout's theorem (cf. Proposition A. 3 of [16j . or [9]), the degree of U is given by the 
formula (|3.4p . which also equals its cardinality. Therefore, the degree of A(/o, /1, . . . , f m ) in 
/o is given by (|3.4p . and then the formula for its total degree immediately follows. 
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c) Clearly, if the singular locus Vp(/i, . . . , f m ) s ing nas positive dimension, then it must 
intersect the hypersurface fo(x) = for arbitrary /o, by Bezout's theorem. Thus the system 
(|3.ip has a singular solution, which implies A(/o, fx, ■ ■ ■ , f m ) = for arbitrary fo- To prove 
the reverse, suppose A(/o,/i, . . . ,/ m ) = is identically zero in /q. We need to show that 
Vp(/i, • • • , f m )sing has positive dimension. For a contradiction, suppose it is zero dimensional 
and consists of finitely many points . . . , u^ N > G P n . Note that the dimension of the set 

T= [J 7J([V 5 /i(u) ••• V £ / m (v)]) 

«6^(/i,...,/m) 

is at most n in the affine space K n+1 whose dimension is n + 1. Here R-(A) denotes the 
column range space of matrix A. So the complement IR n+1 \T has positive dimension, and 
hence we can choose a G M n+1 \T such that the hyperplane a T x = does not pass through 
vfi\ . . . ,u( N \ For f a (x) = a T x, the homogeneous polynomial system 

fa(x) = fl(x) = --- = f m (x) = 

has no singular solution in P", which implies A(/ a , fx, . . . , f m ) ^ by (|3.2p and then 
contradicts that A(/o, fx, ■ ■ ■ , f m ) = is identically zero in /q. So, the singular locus of 
Vp(fx, ■ ■ ■ , fm) must have positive dimension. □ 

The discriminant A(/o, . . . , f m ) of m+1 forms fo(x), . . . , f m {x) is a natural generalization 
of the standard discriminant of a single form. In formula (|3.5p . if we set m = 0, then the 
degree of A(/o) is (n + l)(c?o — l) n , which is precisely the degree of discriminants of forms of 
degree do in n + 1 variables. 

In Theorem 13.11 if every di = d, the discriminant A(/o, . . . , f m ) is homogeneous in every 
fi of degree { m \\)d m (d - I)", and its total degree is (n + l)(^)d m (d - l) n ~ m . This is 
precisely the degree of the discriminant of the multi-homogeneous form C(x, A) (see Theorem 
2.4 of Section 2B in Chapter 13 of 0). 

In Theorem 13. 11 when m = n, the Jacobian of (|3.ip must be singular at its every solution 
u G P™, because by Euler's formula (|2.ip 

^ T [Vx/o(«) • • • V £ f n (u)] = [do/o(fi) • • • dn/»(«)] = 0. 

So (|3.ip has a singular solution if and only if the homogeneous polynomial system 

/o(x) = --- = /„(5)=0 

has a solution in P n , which is equivalent to that the resultant Res(fo, . . . , f n ) vanishes. So 

A(/ ,. ..,/„) = Res(f ,...,f n ) = 0. 

Observe that A(/o, . . . , / n ) and Res(fo, . . . , / n ) have the same degree 

do---d n (dp 1 H h ci" 1 ) . 

So A(/o, . . . , f n ) is equal to Res(fo, . . . , / n ) up to a constant factor. 

When do > 1 and every = /fx (1 < i < m) is linear, (|3.ip has a singular solution if 
and only if fo(x) has a nonzero critical point in the orthogonal complement of the subspace 
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span{/i, . . . , f m }. If every fi{x) = Xi-i, A(/ , x , . . . , x m _i) vanishes if and only if A(/) = 0. 
Here / = /(0, . . . , 0, x m , . . . , x n ) is a form in (x m , ... , x n ). Since A(/ , x , . . . , sc m _i) has 
degree (n — m + l)(do — l) n_ ' m in / , we have 

A(/ ,x ,...,x m _i) = 77- A(/) (3.8) 

for some scalar r\ ^ 0. Furthermore, if /o = x T ^4x is quadratic, then it holds that 

A(x T Ax, xq, . . . , x m _i) = • det A(m + 1 : n + 1, m + 1 : n + 1). (3.9) 

Here -A(7, /) denotes the submatrix of ^4 whose row and column indices are from /. 

We conclude this section by generalizing A(/o, . . . , f m ) to nonhomogeneous polynomials. 
If /o, . . . , f m are not forms, denote by f^ the homogenization of Then A(/o, . . . , f m ) is 
defined to be A(/ ft , . . . , f m ). 

4 Polynomials nonnegative on M. n 

This section studies the cone Pd(K) when K = W 1 . Note that a polynomial f{x) is nonneg- 
ative in R n if and only if its homogenization f h (x) is nonnegative everywhere. So we just 
consider the cone of nonnegative forms. 

Let P n ,d be the cone of forms nonnegative in W 1 of degree d. Here d > is even. Clearly, 
a form / lies in the interior of P n ^d if and only if it is positive definite, that is, f(x) > for 
every i/0. If f(x) lies on the boundary dP n ^, then it vanishes at some 0/m£ W 1 . Since 
f(x) is nonnegative everywhere, u must be a minimizer of f(x) and V/(u) = 0. This implies 
that f(x) has a nonzero critical point, and hence its discriminant A(/) = 0. So the boundary 
dP n ^ lies on the discriminant al hypersurface 

£n,d = {f€ R[x] d : A(/) = 0}. 

Theorem 4.1. The Zariski closure of the boundary dP Ut d is S n ,d> which is an irreducible 
hpyersurface of degree n(d— l) n_1 . 

Proof. The discriminant A(/) is irreducible and has degree n(d — l) n , so the hypersurface 
£ n ^d is also irreducible and has degree n(d — l)™" 1 . Since dP n) d C £ n ,di its Zariski closure 
Zar(dP n) d) lies on £ nt d- The irreducibility of £ n ^d implies Zar{dP n ^d) = £ n ,d- D 

When d = 2, P n ^ reduces to the cone of positive semidefinite matrices. A typical barrier 
for P n ^ is — log det A, where f{x) = x T Ax. Does there exist a similar barrier for P n d when 
d > 2? Unfortunately, this is impossible if we require the barrier to be of log-polynomial 
type, as will be shown in the below. 

Let X m i n (f) denote the smallest value of a form f(x) on the unit sphere 

Knin(f) '■= min f(x). (4.1) 

||!t||2=l 

The boundary dP n ,d is then characterized by \ m in(f) = 0. Clearly, if X m in{f) — then 
A(/) = 0, but the reverse might not be true. For instance, for the positive definite form 
f{x) = \\x\\2 (for even d > 2), \ m i n {f) = 1 but A(/) = 0, because V/(x) = has a nonzero 
complex solution. So the discriminantal hypersurface A(/) = intersects the interior of P n ^. 
when d > 2 is even. This interesting fact leads to the following theorem. 
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Theorem 4.2. If d > 2 is even and n > 2, there is no polynomial </?(/) satisfying 

• tp{f) > whenever f lies in the interior of P n ^d, and 

• <p(f) = whenever f lies on the boundary of P n ,d- 

Therefore, — log </?(/) can not be a barrier function for the cone P U) d when we require (p(f) 
to be a polynomial, and P n d "is not representable by a linear matrix inequality (LMI), that is, 
there is no symmetric matrix pencil 

L(f) = ( Wkere = E ^ 

a£N n :\a\=d a 

such that P n>d = {fe R[x] d : L{f ) y 0} and L(f) y for f G int{P n4 ). 

Proof. For the first part, we prove by contradiction. Suppose such a ip exists. The zero set 
^min(f) = lies on the variety V(ip). Since the discriminant al hypersurface A(/) = is the 
Zariski closure of X m in(f) = 0, i.e., the smallest variety containing \ m in(f) = 0, A(/) = 
is a subvariety of V(<p). So </?(/) is vanishing on A(/) = 0. By Hilbert Nullstenllensatz (see 
Theore UTTj) . there exist an integer k > and a polynomial p(f) satisfying 

p(f) k = A(f)-p(f). 

Now we choose f(x) = \\x\\2 G int{P n ^) in the above, then A(j) = and <p(f) = 0, which 
contradicts the first item. 

For the second part, the non-existence of — log-polynomial type barrier function imme- 
diately follows the first part of the theorem. The non-existence of LMI representation also 
clearly follows the first part, because otherwise the determinant detL(/) would be a polyno- 
mial satisfying the first part. □ 

Theorem 14.21 tells us that there does not exist a polynomial <p(f) such that — log (£>(/) 
is a barrier for P n ,d- However, — log <p{f) would be a barrier if <p(f) is not required to be a 
polynomial. Actually 

</>(/) = " log Xrnin(f) (4.2) 

is a barrier for P nt d, where \ m in{f) is defined by (|4.ip . The function X m i n (f) is semialgebraic, 
positive in int(P nt d), and zero on dP n , t d- The barrier </>(/) is also convex in int{P n ^d)- 

Theorem 4.3. The function 4>{f) is convex in int(P nt d). 

Proof. For any f^\f^' G int(P n! d), from (|4.ip we have 

Xmin + (i - #)/ (2) ) > e\ min (/«) + (1 - 9)\ mm (/w) , V0 € [0, i]. 

Since — log(-) is concave, the above then implies 

4> (ef {1) + (i - #)/ (2) ) < H + (i - (/ (2) ) • 

So 4>(f) is convex in int(P n: d). □ 

However, the barrier — logA m j n (/) is not very useful in practice, because computing 
Amin.(/) is quite difficult. When d = 4, it is NP-hard to compute \ m in(f)- 
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4.1 Computing the discriminantal variety A(/) = 

We have seen that dP n ^ lies on the discriminantal hypersurface A(/) = 0. Cayley's method 
would be applied to compute A(/), as introduced in Chap. 2 of [6]. When n = 2 and d = 4, 
the boundary of i-2,4 lies on the hypersurface defined by the polynomial 

b 2 c 2 d 2 - Aac 3 d 2 - 46 3 d 3 + 18abcd 3 - 27a 2 d 4 - Ab 2 c 3 e + 16ac 4 e + 186 3 cde - 80abc 2 de 
-6ab 2 d 2 e + U4a 2 cd 2 e - 276 4 e 2 + 144a6 2 ce 2 - 128a 2 c 2 e 2 - 192a 2 bde 2 + 256a 3 e 3 , 

where a,b,c,d,e are the coefficients of f(x) = ax\ + bx\x 2 + cx 2 x 2 + dx\x\ + ex\. It is 
a homogenous polynomial of degree 6 in 5 variables. When n = 3 and d = 3, A(/) is 
a homogeneous polynomial of degree 12 in 20 variables, and has 21,894 terms in its full 
expansion. When n = 3 and d = 4, A(/) is a form of degree 27 in 15 variables and has 
thousands of terms. A very nice method for computing discriminants of trivariate quartic 
forms is described in Section 6 of [26] . 

Generally, it is quite complicated to compute A(/) directly. A more practical approach 
for finding the discriminantal locus A(/) = is to apply elimination theory (see [3]). Let 
f p (x) be a form in x whose coefficients are polynomial in a parameter p = (a, b, ...) over the 
rational field Q, i.e., in the ring Q[p]. First, we dehomogenize f p (x) like 

g(l,x 2 , ...,x n ) = f p {l,x 2 , ■ ■ .,x n ). 

If f p (x) £ dP nt d has no nontrivial critical point on the hyperplane x\ = at infinity, then 
the overdetermined polynomial system 

must have a solution. Hence, we can use the elimination method described in [3] to find the 
polynomial equation that the parameter p satisfies. By eliminating x 2 , ■ ■ ■ ,x n m (|4.3p . we can 
get a polynomial tp such that if (|4.3p has a solution then <p(p) = 0. Hence, the discriminantal 
locus A(/ p ) = lies on <p(p) = 0. The polynomial ip(p) = can be found by using function 
elim in software Singular [5] . 

Example 4.4. (i) Consider the polynomials parameterized as 

f a ,b(x) = x\ + x\ + x\ - a{xix\ + x 2 x\ + x$c\) - b(x\x 2 + X2X3 + x\x\). 

Its discriminant ip(a,b) = A(/ a ft) is 

16384(a + 6 - 1) • (a + 6 + 2) 3 • (7a 2 + 76 2 - 13a& + 4a + 46 + 16) 4 - 
(7a 5 + 86a 4 - 17a 4 - 146a 3 + 16a 3 6 2 + 16a 3 - 16a 2 + 486a 2 - 21a 2 6 2 + 16a 2 6 3 
+48a6 2 - 32a6 - 14a6 3 + 8a6 4 - 64a + 76 5 - 176 4 - 166 2 + 166 3 - 646 + 128) 3 . 

The above formula is obtained by using a Maple code that was kindly sent to the author by 
Bernd Sturmfels for computing (3, 3, 3)-resultants. Let 

F = {(a, 6) G R 2 : f afi is SOS in x} . 

It is a convex region in M?. The shape of F would be found by running the following Matlab 
code supported by software YALMIP [13] 
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sdpvar x_l x_2 x_3 a b; 

p = x_l~4+x_2~4+x_3~4-a*(x_l*x_2~3+x_2*x_3~3+x_3*x_l~3) . . . 
-b*(x_l~3*x_2+x_2~3*x_3+x_3~3*x_l) ; 
v = monolist ( [x_l x_2 x_3],2); 
M = sdpvar (length(v) ) ; 

L = [coeff icients(p-v'*M*v, [x_l x_2 x_3] )==0,M>=0] ; 
w = plot(L, [a,b] , [1,1,1] , 100) ; 
fill(w(l, :),w(2, :),'b'); 

The set F is drawn in the shaded area of the upper left picture in Figure [TJ The curves 
there are defined by <p(a, 6) = 0. Since every nonnegative trivariate quartic form is SOS (see 
Reznick [25]), we know F = {(a, b) : f a ,b( x ) £ ^3,4}- 

(ii) Consider the polynomials parameterized as 

/, f ™\ — ™4 ~,4 _i ,~4 _i _i r . j ~, 2 ~, 2 i_ ^2^2 ~, 2 ~, 2 ~,2r„2\ 
a,b\^ J — ~r ^2 ' 3 ' 4 ' ^V^l 2 ' 2 3 x^x^ ^3*^4/ 

+6(x 2 x 2 — x\x\ + X1X2X3X4). 

Eliminating £2,^3; £4 in (|4.3|) gives 92(0,6) as 

(a + 2) • (a - 2) • (b + 2) • (b - 2) • (16a 2 + 16a6 + 5b 2 + 32a + 166 + 16)- 
(16a 2 - 16a& + 5& 2 - 32a + 166 + 16) • (4a 2 6 - 8a 2 - 56 2 + 16)(56 2 - 166 + 16). 

The curve A(/ a) b) = lies on ip(a, 6) = 0. Let 

F = {(a, 6) G R 2 : f afi is SOS in x} . 

It is a convex region. Using the method in (i), we get F is the shaded area of the upper right 
picture in FigureHJ The curves there are defined by ip(a, 6) = 0. Let G = {(a, 6) : / 0) & G Pi,4}- 
Clearly, F C G and the boundary of G lies on ip(a, 6) = 0. From the picture, we can see that 
F is a maximal convex region whose boundary lies on (p(a, 6) = 0. So F = G. 

(iii) Consider the polynomials parameterized as 

fa,b{x) = x\ + x\ + x% - a{x\(x\ + xf) + x\{x\ + x\) + x\(x\ + x^fj + bx\x\x\. 

When a = 1, 6 = 3, /i,3(x) becomes Robinson's polynomial that is nonnegative but not SOS 
(see Reznick |25j). Robinson's polynomial has 10 nontrivial zeros, so /1 3 G -P3,6- Eliminating 
X2,X3 in (|4.3p gives tp(a,b) as 

(a - 1) • (a + 3) • (3a + 6 + 3) • (6a - 6 - 3) • (2a 3 + a 2 6 + 3a 2 - 6 2 + 36 - 9). 

The curve A(f a ^) = lies on ip(a, 6) = 0. Let 

F = { (a, 6) G R 2 : {x\ + x\ + x\ + a£)/o,6 is SOS in x) . 

It is an unbounded convex set in R 2 . To get the shape of F, we bound a, 6 as a + 5 > 0, 
40 — 6 > 0. Using the method in (i), we get F is the shaded area of the lower left picture in 
Figure [U The curves there are defined by ip(a, 6) = 0. Let G = {(a, 6) : f a ^ G P^fi}- Clearly, 
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Figure 1: The pictures of curves (p(a,b) = and regions F for polynomials f a ,b( x ) in 
Example 14,41 The upper left is for (i), the upper right for (ii), the lower left for (iii), and the 
lower right for (iv). 

F C G and the boundary of G lies on ip(a, 6) = 0. If f a ,b(x) £ P3fi, then / a ,&(l, 1, 1) > and 
/^(l, 1,0) >0 imply 

b > 6a — 3, a < 1. 

From the picture, we can see that F is a maximal convex region whose boundary lies on 
(p(a, b) = and satisfies the above two linear constraints. So F = G. 
(iv) Consider the polynomials parameterized as 

/i ( <~r*\ — ( t»2 i i ry 2 \ 2 ( ™2™,2 i ™2™,2 i ^y.2™,2 i ™2™,2 i ™2™,2\ 

a,6l x ; — l x l "I" ~r x 5^ 0.^x^X2 "I" X 2 X 3 -|- X3X4 -r X 4 X 5 -f x 5 x ij 

-6(xf + x| + x\ + x\ + xf ). 

When a = 4, 6 = 0, f^o(x) becomes Horn's polynomial (see Reznick |25j). Eliminating 
X2,X3,X4,X5 in (I4.3P gives tp(a,b) as 

(a + 6 - 5) • (a - 2b) ■ (a + 26 - 4) • (b - 1) ■ b ■ (b - 2) • (a 2 + 2ab - 46 2 )- 
(a 2 - 26 2 - 4a + 66) • (a 2 - 2a6 - 46 2 - 4a + 166) • (a6 + 26 2 - a - 66). 

The curve A(/ Q) &) = lies on 99(0, 6) = 0. Let 

F = {(a, 6) G R 2 : (x 2 + x 2 . + x\ + xf, + xD/^ft is SOS in x} . 

It is also an unbounded convex set. To get the shape of F, we bound a,6asa+2 > 0, 6+4 > 0. 
Using the method in (i), we get F is the shaded area of the lower right picture in Figured! 
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The curves there are defined by (p(a,b) = 0. Let G = {(a, b) : f a ,b G -P3,6}- Clearly, F C G 
and the boundary of G lies on <p(a,b) = 0. Then / 0>6 (1, 0, 0, 0, 0) > 0, / 0)6 (1, 1, 0, 0, 0) > 0, 
/ aj 6(l, 1,1,1,1) > imply that any pair (a, b) £ G satisfies 

a + 6-5<0, a + 26-4<0, 6 - 1 < 0. 

Since (3.10, 0.5), (5.5, -1), (9.1, -4) (£ G (verified by software GloptiPoly 3 [IT), by observing 
the lower right picture in Figure [TJ we can see that F is a maximal convex region that satisfies 
the above three linear constraints, excludes the previous 3 pairs, and has the boundary lying 
on (p(a,b) = 0. So F = G. □ 

4.2 Nonnegative multihomogeneous forms 

In this subsection, we study the cone of nonnegative multihomogeneous forms. Let M^''"'^ r 
denote the space of multihomogeneous forms in the space M. ni x • • • x M nr which are homo- 
geneous of degree di in each M. di . Thus every / € M^'"''^ r has the form 

/= E / ai a.f^r-^r. 

(ai,...,a r )€N"ix-xN"'- 

Here we assume all the degrees di are even. Let P^ 1 ''"'^ be the cone of forms in M^'"'^ r 
that are nonnegative everywhere. 

Given / e M^ 1 ;//^ , we say (uW , . . . , «( r )) € ]J r i=1 C" ! is a critical point of / in ]J r i=1 P"*" 1 
if every ^ and 

V x(1) f(u ( - 1 \...,u^) = 0, V x(r) /(^ 1 ...,nM) = 0. 

Let //; ; ::;;; C be the set 

H Z::S: = {/ ; M Z::t ■ f has a critical p° int in II '}• 

It was shown in [6l Prop. 2.3 in Chap. 13] that H^''"'J is a hypersurface if and only if 

2(rii - 1) < ni H h n r - r for all i: di = 1. (4.4) 

In particular, if every di > 1, H^'"''^ r is a hypersurface for any dimensions ri\, . . . , n r . When 
(|4.4p holds, we still denote by A(/) a defining polynomial of the lowest degree for H^''"'^ . It 
can be chosen to have coprime integer coefficients and is unique up to a sign. The polynomial 
A(/) is also called the discriminant of the multihomogeneous form /. 

Theorem 4.5. When all di > are even, the boundary dP^''"'^ lies on the hypersur- 
face H^i ' dr w h° se degree is the coefficient of the term • • • z™ r_1 in the power series 
expansion of the following rational function 



djZj 



J =1 
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Proof. Since all di > are even, the condition (|4.4p holds, and H^''"'^ is a hypersurface 
defined by A(/) = 0. A multihomogeneous form / E P^'" if and only if 

A mm (/):= min /(x«,..,x«) > 0. 

||x( 1 )|| 2 =--'=||2;M||2=l 

Clearly, / G dP^^Z if and only if A ™«(/) = °- If f G dP Z''Sr > then we can find 
uW, . . . , tt^ of unit length satisfying f(yP~\ . . . , u^) = and 

V a( i)/(«W,...,«W) = 0, V x(r) /(«W J ...,«W) = 0. 

Thus, / also belongs to H^}'"'2 r . The degree formula for H^''"'^ is given by Theorem 2.4 
of Chapter 13 in [6]. □ 




Example 4.6. Consider the bi-quadratic forms parameterized as 

fa,b{x) = {x\ +x\ + xl){x\ + x\ + xl) + a(x\x\ + x\x\ + x\xl) 

+&(xiX 2 X 4 X 5 + £1X3X4X6 + X 2 X3X 5 X 6 ). 

Here ni = 712 = 3, d± = d2 = 2. First, we dehomogenize f a ,b{x) as 5 = / a ,&(l, X2, X3, X4, 1, Xe), 
and then use the function e/im in Singular to determine all pairs (a, 6) satisfying A(f a ^) = 0. 
Eliminating X2,X3,X4,X6 from 

<9g <9g dg dg 

® 8x2 8x3 8x4 Oxq 

gives the equation ip(a, b) = where p(a, b) is 

(a + 1) • (a + 6 + 3) • (a 2 - ab + 6 2 ) • (-6 2 + 4a + 46) • (-0 2 + 4a - 46)- 

(a 3 o 4 - 16a 6 - 8a 4 o 2 - 4a 3 6 3 + 3a 2 6 4 + ab 5 - 80a 5 - 16a 4 6 - 32a 3 6 2 - 20a 2 6 3 

-4ao 4 + b 5 - 96a 4 - 32a 3 6 - 24a 2 6 2 - 12a6 3 - 5o 4 ). 
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The curve A(f a ^) = lies on ip(a, b) = 0. Let 

F= {(a,b) G R 2 : (1 + x\ + x\ + x\ + xf) ■ / a , b (l, x 2 , x 3 , x 4 , 1, x 6 ) is SOS } . 

By the method used in Example 14. 44 F is drawn in the shaded area of Figure [2j The curves 
there are defined by (f(a,b) = 0. Let G = : f a ,b(x) G -P^}' Clearly, F C G and the 

boundary of G lies on <^(a, 6) = 0. If f a bi x ) £ ^2 2' then from 

/ a , 6 (l, 1, 1, 1, 1, 1) > 0, / a , 6 (l, 0, 0, 1, 0, 0) > 

we know every (a, b) £ G satisfies 

a + 6 + 3 > 0, a + 1 > 0. 

Because /2o,i5(^) G (•.' V^ 1)X2)a . 3 /2o,i5 has negative eigenvalue at (1, 1,0)) and /2o,-is(^) 
G ('■' x 2 3:3/20,-15 has negative eigenvalue at (1, —1,0)), from the picture we can see that 
F is a maximal convex region that excludes (20, 15) and (20, —15), satisfies the above two 
linear constraints, and has boundary lying on ip(a, b) = 0. So F = G. □ 

5 Polynomials nonnegative on a variety 

This section studies the cone Pd(K) when K is a real algebraic variety defined as 

K = {x G R n : gi(x) = ■■■= g m (x) = 0}. 

Here g = (gi, . . . ,g m ) is a tuple of polynomials. For convenience, denote Pd{K) as 

Pd{g) = {fix) G R[x]< d : f{x) > for every x G V R (g)}. 

To study the boundary dPdig) of Pdig), we need a characterization for it. One would 
think if / lies on dPdig) then fix) vanishes somewhere on Vk(#). However, this is not always 
true. For a counterexample, consider / = x\ + x 2 and g = x\ + x 2 — 1- Clearly, / is strictly 
positive on Mr(#), but it lies on dP±ig). For any e > the polynomial x\ + x 2 — e is no longer 
nonnegative on Mr(<?) because 

inf x\+x 2 = 0. 

The reason is that V^ig) is not compact. We need other characterization in this case. 
Let V^ig) be the homogenization of V^ig), that is, 

V£{g) = [x G R n + i . = . . . = = o j . 

Clearly, if / ft is nonnegative on V£(g), then / is also nonnegative on Vr (g), but the reverse 
is not necessarily true. For this purpose, we need a new condition. We say the variety Vj^(g) 
is closed at 00 if 

V£ig) n {x > 0} = closure (y£(g) D {x > 0}) . 
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Define two constants 

S g (f):= min f(x), (5.1) 
xevk(g) 

8 h g {f):= min f\x). (5.2) 

xey^(g):||x||2=l,xo>0 

The boundary dPd(g) is characterized as below. 
Proposition 5.1. Let g be given as above, 
(i) IfV$i{g) is compact, then 

S g (f)>0 & f €int(P d (gj), and S g (f) = O / G 9P d ( 5 ). 

(m,) IfV^(g) is closed at oo, then 

5 h g (f)>0 & f €int(P d (g)), and S^(f) = ^ f e dP d (g). 

Proof. Part (i) is quite clear. We prove part (ii). For any u G V^(g) with uq > 0, we can 
find a sequence (th,w k ) G V^(g) with every ^ > approaching u. Note that w k /t k G Vr (g). 
So, if / G Pd(p), then 

f h (u) = lim / h (t fc ,^ fc ) = lim 4f(w k /t k ) > 0, 

and we have <^(/) > 0. On the other hand, if Sg(f) > 0, then for every jj G 14(g) 

/(t/) = / h (l,t;) = (1 + \Hlf 2 f h ((l,t/)/(l + ll^li) 1 / 2 ) > (1 + \Hlf 2 S^f) > 0, 

and we get / G P d {g). The above implies 5 g (f) > if and only if / G P d {g). 

By definition, 6g(f) is the minimum of a polynomial function over a compact set. If 
Sg(f) > 0, then in a small neighborhood O of / we have 5 g (p) > for every p G 0, that is, 
/ lies in the interior of P d (g). If $g(f) = 0, then we can find p G M[x]<rf of arbitrarily small 
coefficients such that 5 g (f + p) < 0, that is, / G dP d (g). □ 

We would like to remark that not every V£(g) is closed at oo, and even if Vu.(g) is compact 
Vfil(g) might still not be closed at oo. 

Example 5.2. (i) Let g = x\ {x\ — x%) — 1 and / = x\ — Xi + 1. The polynomial / is strictly 
positive on the variety Vs.(g), but f h = x\ — X2 + xq is not nonnegative on 

Vr(9) = {(x ,x 1 ,x 2 ) G M 3 : x\{x 1 -x 2 )-xl = 0} . 

This is because (0,0, 1) G V^(g) while f h (0,0, 1) < 0. So V£(g) is not closed at oo. 
(ii) Let g = x\ (1 — x\ — x 2 ) — x\. The variety V^{g) is compact. Its homogenization is 

V£(g) = {x : xj(x 2 -x\- x\) - x\x\ = 0} . 

However, V$(g) is not closed at oo. Otherwise, for every u G V£(g) n {xq = 0} we have 

u= lim t k (l,v k ) for some u fe G Vi(st). 

tfc>0,t fc ->0 

This implies V£(g) n {xq = 0} is compact, which is clearly false. □ 
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Now we study the boundary of the cone Pd{g)- 
Theorem 5.3. Let g = (gi, . . . ,g m ) be given as above, and deg(^) = d{. Suppose m <n. 

(i) If Vr(<7) 0, and either Vr(<?) is compact or V^(g) is closed at oo, then the boundary 
dPd(g) lies on the hypersurface 

£ d {g) = {/ e R[x]< d : A(f, 9l , ...,g m ) = 0}. 

(ii) If the projective variety Vp(gi, . . . , <? m ) ^ s nonsingular, the degree of £d(g) is 

f n d * \ ■ Sn ~ m ( d ~ ij d ~ i,di ~ ^ • • • ' dm ~ -0 • < - 5,3 ) 

Otherwise, the above is only an upper bound. 

(Hi) The polynomial A(f,gi,...,g m ) is identically zero in f if and only if the projective 
variety Vp(gY , . . . , g m ) has a positive dimensional singular locus. 

Proof, (i) We first consider the case that V^(g) is closed at oo. Let f{x) £ dPd(g). By 
Proposition 15.11 we know f h is nonnegative on V^(g) and vanishes at some 0/ii£ V^(g). 
So u is a minimizer of f h (x) on By Fritz-John optimality condition (see Sec. 3.3.5 in 

[T]), there exists (/j,q, m, . . . , /i m ) 7^ satisfying 

[J>oVxfo(u) + mVxgi(u) H h HmV^gmiu) = 0, 

/C") = 51 («) = ••• = ffm(w) = 0. 

By relation (|3.2p . we know A(f,g±, . . . ,g rn ) = 0. 

The proof for the case that Mr (5) is compact is almost the same as the above, and is 
omitted here. 

(ii) When Vf{g\, . . . , <7 m ) is nonsingular, from the proof of part b) in Theorem 13. 1\ we 
know the degree of A(/, gi, . . . ,g m ) in / is given by (15. 3j) . When Vp(gi, ■ ■ ■ , <? m ) i s singular, 
the formula in (|5.3p is only an upper bound by perturbing the coefficients of g±, . . . ,g m . 

(iii) This immediately follows part c) of Theorem 13.11 □ 

We have seen that there is no log-polynomial type barrier function for the cone P^(M n ) 
when d > 2 and n > 1. There is a similar result for Pd(g). 

Theorem 5.4. Suppose V5r((?) is nonempty, either V^(g) is compact or V^(g) is closed at 
00, Vf(g h ) has positive dimension, and d > 2 is even. If the discriminant A(/, g±, . . . ,g m ) is 
irreducible in f over C, then there is no polynomial </?(/) satisfying 

• (f(f) > whenever f lies in the interior of Pd(g), and 

• (f(f) = whenever f lies on the boundary of Pd(g). 

Therefore, —log(p(f) can not be a barrier function for the cone Pd{g) when we require </?(/) 
to be a polynomial, and Pd(g) is not representable by LMI. 



19 



Proof. We prove the first part by contradiction. Suppose such a (p exists. By Theorem 15.31 
we know dPd(g) lies on the hypersurface A(f,gi,...,g m ) = 0. Since A(/, g±, . . . , g m ) is 
irreducible in /, the hypersurface A(/, gi, . . . ,g m ) = is irreducible and equals the Zariski 
closure of dPd(g) (it is contained in some hypersurface). Hence, the hypersurface tp{f) = 
contains A(/, g±, . . . , g m ) = 0, and (f(f) vanishes whenever A(/, g±, . . . , g m ) = 0. By Hilbert's 
Nullstenllensatz (see Theorem 12. lj) . there exist an integer k > and a polynomial p(f) such 
that 

<p{f) k = A(f, gi ,...,g m ).p(f). 

Set f{x) = (1 + x\ -\ h x 2 n ) d l 2 , then f h (x) = (x% + x\ -i h x 2 n ) d / 2 . Clearly, / lies in the 

interior of Pd{g)- However, since Vp(g h ) has positive dimension, we know 

Vw(f h ,g h ) = {xeV n :xl + x 2 1 + --- + x 2 n = 0}nV(g h )^® 

by Bezout's theorem. For any u G Vp(f h ,g ), we have V^/ /l (n) = (d > 2) which results in 
A(/, gi, . . . ,g m ) = 0. So (p{f) = 0, which contradicts the first item. 

The second part is a consequence of the first part, as in the proof of Theorem 14.21 □ 

5.1 Computing the discriminantal variety A(f,gi, . . . , g m ) = 

Now we discuss the connection between A(f,gi,...,g m ) and the discriminant of the La- 
grangian polynomial in (x, A) 

k 

L{x, A) = f(x) + ^2 ^i9i{%)- 

i=l 

When Vr(<?) is compact, / E dPd(g) if and only if 5 g (f) = 0, i.e., there exists u G Vr(s() 
such that f{u) =0 and u is a minimizer of / on Vu.(g). So, if f{x) G dPd(g) and Vu.(g) is 
nonsingular at u, the Karush-Kuhn- Tucker (KKT) condition (see Sec. 3.3 in [lj)holds, and 
there exists /i = (/ii, . . . , fi m ) satisfying 

m 

Vx/(u) + y~jtijV x9i{u) = 0, (u) = • • • = g m (u) = 0. 

i=l 

The above is equivalent to that (u, /i) is a critical zero point of L(x, A), that is, 

V X) Ai(u, M) = 0, L(u, /i) = 0. 

Hence, we have A(L) = 0. Therefore, the hypersurface A(/, g\, . . . , g m ) = would be 
possibly determined via investigating A(L) = 0. To the best knowledge of the author, no 
general procedure is known in computing the discriminant of type A(f,g±, . . . ,g m ). Though 
there exist systemic methods for evaluating A(L), its computation and formula would be too 
complicated to be practical, as we have seen in the preceding section. In the following, we 
propose a different approach using elimination. 

Suppose / = f(x;p) is a polynomial in x whose coefficients are also polynomial in a pa- 
rameter p = (a, 6, c, . . .) over the rational field, i.e., from the ring Q\p\. So, if f(x;p) G dPd(g) 
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and Vr (g) is a nonsingular compact set, then / satisfies the over-determined polynomial 
system in (x, A) 

m \ 

V*/(z) + Y,\iV x 9i(x) = I 

/(z) = = • • • = 5m(a;) = J 

The equation that p satisfies would be determined by eliminating (x, A) in the above. Let 
<p(p) = be the polynomial equation obtained by eliminating (x, A) in (15. 4p . So, if p satisfies 
A(/, gi, . . . , g m ) = 0, then cp(p) = 0. Computing <^(p) would be done by using elim in 
Singular [8]. We illustrate this in the below. 

Example 5.5. Consider the polynomials parameterized as 

/ = x 2 + ax\X2 + bx\ + CX2 + d, 

and -fT = {x 2 + x?> = 1} is a circle. The polynomial <£>(a, 6, c) obtained by eliminating (x, A) 
in ([53D is 

a 6 _ 3a 4 6 2 + 3a 2 6 4 _ 6 6 _ 3a 4 c 2 _ 2 l a 2 b 2 c 2 - 36 4 c 2 + 3a 2 c 4 - 36 2 c 4 - c 6 
+36o 3 6cd + 18ab 3 cd + I8abc 3 d - 8a A d 2 - 20a 2 b 2 d 2 + b 4 d 2 - 20a 2 c 2 d 2 
+2b 2 c 2 d 2 + c 4 d 2 - Wabcd 3 + 16a 2 d 4 + 18a 3 bc - 18ab 3 c + 36a6c 3 - 8a 4 d 
-2a 2 b 2 d + 106 4 d - 38a 2 c 2 d + 26 2 c 2 (i - 8c 4 d - 24abcd 2 + 32a 2 d 3 - 8b 2 d 3 
+8c 2 d 3 + a 4 - 2a 2 6 2 + b 4 - 20a 2 c 2 + 206 2 c 2 - 8c 4 + 2Aabcd + 8a 2 a? 2 - 326 2 d 2 
-8c 2 d 2 + 16(i 4 + 16a6c - 8a 2 d - 8b 2 d - 32c 2 d + 32d 3 - 16c 2 + 16d 2 . 

It is a polynomial of degree 6 in 4 variables. The set {(a,b,c) : / G clP2(i^)} lies on the 
surface (p(a, 6, c) = 0. □ 

Example 5.6. (i) Consider the polynomials parameterized as 

/ = x 4 + axfx2 + bx\x\ + c, 

and .FT = {xf + x\ = 1} is a circle. The polynomial (p(a,b,c) obtained by eliminating (x, A) 
in (pT4"l) is 

4a 3 6 3 + 27a 4 c 2 - 36a 3 6c 2 + 2a 2 6 2 c 2 - 36a6 3 c 2 + 276 4 c 2 - 256a 2 c 4 
+512a6c 4 - 25 66 2 c 4 + 6a 2 6 2 c - 36a6 3 c + 546 4 c - 288a 2 c 3 + 704a6c 3 
-5446 2 c 3 + 276 4 + 192a6c 2 - 2886 2 c 2 - 256c 4 - 256c 3 . 

The surface (f(a, b,c) = is drawn in the left picture in Figure [3l It contains the set 
{(a,b,c) :fedP 4 (K)}. 

(ii) Consider the polynomials parameterized as 

/ = x\ + ax 3 X2 + 6x1X2 + c, 

and K = {xf + x\ = 1} is a circle defined in 4-norm. The polynomial <p(a, b, c) obtained by 
eliminating (x, A) in (|5.4p is 

4a 3 6 3 + 27a 4 c 2 + 6a 2 6 2 c 2 + 276 4 c 2 + 192a6c 4 - 256c 6 + 6a 2 6 2 c 
+546 4 c + 384a6c 3 - 768c 5 + 276 4 + 192a6c 2 - 768c 4 - 256c 3 . 

The surface (p(a, b,c) =0 is drawn in the right picture in Figure [3l It contains the set 
{(a,b,c) : f eOP 4 (K)}. 

The surfaces in Figure [3] are drawn by Labs' software Surfex which is downloaded from 



the website |www . surfex . algebraicsurf ace . net □ 
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Figure 3: The pictures of surfaces if (a, b,c) = in Example 15.61 The left is for (i), and the 
right is for (ii). 

5.2 Resolution of singularities 

In Theorem 1 5, 3 1 we know if the projective variety Vp(g h ) has a positive dimensional singular 
locus, then A(/, gi, . . . , g m ) is identically zero in / and A(/, g\, . . . , g m ) = defines the whole 
space M. n . This is not what we want, because the boundary dPa(g) usually has codimension 
one. To study dPd(g), we need to resolve the singularities of Vp(g h ). By Hironaka's result 
(see Theorem 17.23 in Harris's book [9]), there exist a smooth projective variety U C P n and 
a rational mapping 

<j> : U — > V F (g h ) 

such that <p(U) is dense in Vp(g h ). Thus, / G Pd{g) if an d only if f h (<fi) is nonnegative on U. 
Consequently, the boundary of Pd(g) can be investigated through studying forms nonnegative 
on U. We illustrate how to do this as below. 

Example 5.7. Consider the variety V(g) C C 3 where 

g(x) = (( a;i -l)2 + x 2_ l) 3_ x | 

Both V(g) and V^(g h ) have positive dimensional singular locus. Let 

U = {y G P 3 : y\ + y\ - y yi - yt = 0}. 
It is a smooth variety. Let <f> be the mapping: 

4>- y = (2/0,2/1,2/2,2/3) 1 — > x = (2/0,2/0 + 2/1, 2/2,2/3)- 

Then <j>{U) = V P (g h ). So, f(x) G P d (g) if and only if f\<j>) G P 3d (g), and /(x) G 9P d (<7) if 
and only if f h ((j)) G dP 3d (q). Here g = yf + y\ - y 2/| - 2/o- 1=1 

However, we would like to remark that such eft and U are typically quite difficult to find. 
This issue is beyond the scope of this paper. 



22 



6 Polynomials nonnegative on a semialgebraic set 

This section studies the cone P d (K) when K is a general semialgebraic set in W 1 . Consider 
K is given as 

K = {x£R n : ( 5 i (x), . . . ,g m (x)) = 0, ( Pl (x), . . . ,p t (x)) > 0}. 

Here the gi and pj are all polynomials in x. Recall that 

P d {K) = {/ G R[x]< d : f(x) > Vx G K}. 

We are interested in the algebraic geometric properties of its boundary dPd(K). Typically, 
it is a union of hyper surf aces. 

We begin with the characterization of the boundary dP d (K). Like the case of K being a 
real algebraic variety, a polynomial positive on K may not lie in the interior of P d (K). Let 
K h be the projectivization of K which is defined as 

K h = [x € M n+1 : • • • , (£)) = 0, (pf (£), . . . > o} . 

Define two constants 

6 K (f) = mm /(s), (6.1) 

5 h K (f) = min f\x). (6.2) 

x€X' l :||x||2=l,a;o>0 

Similarly, we say K h is closed at oo if 

K h n {x > 0} = closure (K h n {x > 0} 

We would like to remark that the definitions of K h and <%(/) depend on the defining poly- 
nomials of K that are usually not unique. So in the places where K h or 5j ( (f) appears, we 
usually assume the defining polynomials of K are clear from the context. 

The interior and boundary of the cone Pd{K) are characterized in the proposition below, 
whose proof is almost the same as for Proposition 15.11 

Proposition 6.1. Let K be given as above. 

(i) If K is compact, then 

M/) > / G int(P d (K)), and 6 K (f) = / G 5P d (K). 

(m,) If K h is closed at oo, i/ien 

>0 ^ f eint(P d (K)), and 5 h K (f) = O / G 5P d (K). 

Using the above characterization, we can get the following result about dP d {K). 



23 



Theorem 6.2. Let K be given as above. Assume at most n — m inequality constraints are 
active at any nonzero point in K h . If either K is compact or K h is closed at oo, then the 
boundary dPd(K) lies on the hyper surface 

£ d (K) := i / G R[x]< d : J] A(f, 9l , . . . ,g m ,p h , . . . ,p ik ) = 

^ {ii,...,ife}C[i],fc<n— m 

Proof. Let f(x) 6 dPd(K). First assume K h is closed at infinity. So there exists 0^m£ K h 
such that f h {u) = 0. Let . . . ,i^} be the index set of active inequality constraints 

= ■■■=!£(«) = o. 

By assumption, k < n — m. Note that u is a minimizer of / fe on K h . By Fritz- John optimality 
condition (see Sec. 3.3.5 in [TJ), there exists Qzo,Mi> • • • , ^ m +k) 7^ satisfying 

m k 

^xf h (u) + EwV^fW + EAWiVij>J(«) = 0, 

/ h (u) = <#(u) = • • • = = p£(u) = • • • = tf h (u) = 0. 

So u is a singular solution to the polynomial system 

f\x) = g1{x) = ■■■= g h m (x) = P l{x) = ...= p^x) = 0. 

Hence, A(f,g 1 ,,...,g m ,p il ,...,p ik ) = 0. 

The proof is similar when K is compact. □ 




-4-3-2-1 1 2 3 4 
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Figure 4: The picture of ip(a, b) = and the set F in Example I6.31 
Example 6.3. Consider the polynomials parameterized as 

fa,b(x) = x\ + x\ + a(xfx 2 + X\x\) + + £2) + 1, 

and K = {1 - x\ - x\ > 0} is a ball. From Theorem 16.21 the boundary of P^{K) lies on the 
union of A(f a ^) = and A(f a) f,,g) = 0. The discriminant q(a,b) = A(/ a> &) is 

2097152(a + l) 2 (a - l) 3 (a 2 + 8) 4 (32 + 32a - 276 4 )(256 + 32a 2 + 276 4 - 27a6 4 ) 2 . 
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By the method used in subsection 15.11 eliminating (x, A) in (|5.4p gives h(a, b) = where 
h(a, b) is 

(a + 2V2b + 3) • (a - 2^2b + 3) • (a 5 + a 3 b 2 - 3a 4 
-30a 2 o 2 - 27b 4 + 32a 3 + 48a& 2 - 96a 2 + 2246 2 + 256a - 768). 

The curve A(/ 0i 6,g) = lies on h(a,b) = 0. Let ip(a,b) = h(a,b) • q(a,b). The curves in 
Figure H] are defined by <p(a,b) = 0. Let 

fa,b( x ) = a o( x ) + °1 0*0(1 - INI!) 
(a, b) E R 2 : cr (x), ^(x) are SOS in x 

deg(ao) = 4, deg(cri) = 2 

It is clearly a convex set. By the method used in Example 14.41 P 18 drawn in the shaded 
area of Figure [H Let G = {(a,b) : f a ^ E Pa(K). Clearly, F C G and the boundary of G lies 
on ip(a,b) = 0. Since the polynomials /2,i.5, /2,-i.5) /4,o are n °t nonnegative on the unit ball 
(verified by GloptiPoly 3 [II]), we know (2, 1.5), (2, -1.5), (4,0) G. From Figure HI we can 
observe that F is a maximal convex region that excludes the pairs (2, 1.5), (2, —1.5), (4, 0) 
and has the boundary lying on (p(a, b) = 0. So F = G. □ 



Now we discuss the barriers for P^K). The following is similar to Theorem 14.21 

Theorem 6.4. If K has nonempty interior, d > 2 is even and n > 1, then there is no 
polynomial <p(f) satisfying 

• (f(f) > whenever f lies in the interior of Pd(K), and 

• <£>(/) = whenever f lies on the boundary of Pd(K). 

So, —log (f(f) can not be a barrier function for the cone P^{K) when we require <£>(/) to be 
polynomial in f , and Pd{K) is not representable by LMI. 

Proof. Prove by contradiction. Suppose such a (p exists. Since int(K) / 0, one piece of the 
boundary dPd(K) must lie on the irreducible discriminantal hypersurface A(/) = 0. The 
rest of the proof is then almost the same as for Theorem 14.21 and is omitted here. □ 

Typically there is no log-polynomial type barrier for the cone Pd(K). However, Pd(K) has 
log-semialgebraic type barriers. When K is compact, — log 6k(/), or when K h is closed at oo, 
— log Sf^(f), is a convex barrier for Pd(K), because both Sxif) and 5^(f) are semialgebraic, 
positive in int{Pd{K)), zero on dPd(K), and concave in /. Generally, it is quite difficult to 
compute 6x(f) or £#(/) for general / and K. So these two barriers are not very useful in 
practice. 



6.1 Co-positive polynomials and matrices 

A form f(x) is said to be co-positive if f(x) > for every x E . Clearly, f{x) is co-positive 
if and only if its associated even form 

q f (x) = f(x1,...,x 2 n ) 

is nonnegative in R n . A symmetric matrix A is called co-positive if the associated quadratic 
form f(x) = x T Ax is co-positive. 
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Let C nt d be the cone of copositive forms in R[x]d, and dC n4 be its boundary. Clearly, if 
/ € dC n) d, then there exists / ti 6 E" such that /(u) = 0, or equivalently qf{y/v) = 0. 
Thus dC n< d lies on the discriminant al hypersurface A(gj) = 0. 

Proposition 6.5. The Zariski closure of dC ni d is the hypersurface 
£ d (Rl) :=lf€ R[x] d : ]J A(f I (x I )) = 

{ ^IC[n] 

Here xj = (xj : i € /) and // is obtained from f(x) by setting Xj = for j /. 

Proof. Let / G dC n ^. Then there exists ^ « £ W] such that /(it) = 0. The index set 
/ = {i : m > 0} C [n] is nonempty, and fi(xj) has a positive critical zero point, because 
//(uj) = 0. So A(//(x/)) = 0. Hence, we have Zar(dC n ^) Q £d(]R™). To prove they 
are equal, we need to show that A(//(x/)) = lies on Zar(C n ,d) for every ^ I C [n]. Fix 
such an arbitrary /. Let fi(xj) be a co-positive form which vanishes at 1/ (1 is the vector 
of all ones). Then there is a neighborhood U of // such that every gj £ U n Cu\d vanishes 
somewhere near 1/. Thus U flCin ^ C {A(fj(x[)) = 0}, and 

Zar{U n C|j| i(i ) C Zar({A(fj(xj)) = 0}) = {A(fj(xj)) = 0}. 

Since the hypersurface A(fj(xj)) = is irreducible, we must have 

{A(/;(x/)) = 0} C Zar(UnC m ) C Zar{C n4 ). 

The above is true for every ^ I C [n]. So Zar{C n4 ) = f^(R"). □ 

Proposition 16.51 is equivalent to the fact that 

A( 9/ ) = A(/ 7 ) = 0. 

This is because V x (Qf(x)) = 2diag(x) • V x f(x 2 ) and 

A(g/) = Res • • • , (z 2 )) = 



n a (^)=°- 

0^/C[n] 

We refer to Theorem 1.2 in [6j Chapt.10] for the last equivalence in the above. If [n]\I = 
{i±, . . . , ifc}, (|3.8p implies A(//(x/)) = r]A(f,Xi 1 , . . . ,Xi k ) for some rj 7^ 0. In particular, if 
d = 2 and /(x) = x T ^4x is quadratic, then Proposition 16.51 and (|3.9p imply Zar(dC n ^) is the 
hypersurface 

]J det A(/, /) = 0. (6.3) 

0^/C[n] 
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Corollary 6.6. Suppose d>2 and n > 2. Then there is no polynomial </?(/) satisfying 

• <p(f) > whenever f is in the interior of C n< i, and 

• ip(f) = whenever f is on the boundary ofC n( ^. 

So, —log <p(f) can not be a barrier function for the cone C n ^d when we require <p(f) to be 
polynomial in f , and C n d is not representable by LMI. 

Proof. Prove the first part by contradiction. Suppose such a </?(/) exists. Then 

¥>(/) = o v/ecKV 

So the Zariski closure of dC n d lies on the hypersurface <p(f) = 0. Since d > 2, A(/) is an 
irreducible polynomial in /. By Proposition 16.51 t ne hypersurface A(/) = lies on </?(/) = 0, 
and <p(f) vanishes on A(/) = 0. By Hilbert Nullstellensatz (see Theoren l2.ip . there exist a 
positive integer k > and a polynomial </>(/) such that 

^(/) fc = 0(/)A(/). 

In particular, if we choose / to be f(x) = (l^x) d in the above, then 

<p(f) k = <f>(f)A(f) = 0. 

This is because the form f(x) has a nonzero critical point when d > 2 and n > 2. However, 
f(x) clearly lies in the interior of C n> d, which contradicts the first item. 

The second part clearly follows the first part. □ 

Remark: Corollary 16.61 would be implied by Theorem 16.41 for the case that d > 2 is even. 

Example 6.7. (i) Consider the symmetric matrices A parameterized as 



A = 



1 


a 


-b 


b 


a 


1 


-b 


—a 


b 


-b 


1 


—a 


b 


—a 


—a 


1 



We are interested in the set of all pairs (a, b) such that A is co-positive. The polynomial 
</j(a, b) defining equation (|6.3p is 

-(a - l) 5 • (a + l) 3 • (6 - l) 3 • (b + l) 5 • (-2b 2 + a + l) 2 • (2a 2 + b - l) 2 - 
(a 2 + 3ab + a + b 2 - b - l) • (-a 2 + ab + a - b 2 - b + l) . 

The curve tp(a,b) = is drawn in the left picture of Figure [5j Let 

F = {(a, b) £ R 2 : A = X + Y,X t 0,Y > 0} . 

By the method used in Example 14. 4\ F is drawn in the shaded area of the left picture in 
Figure [5j Because every co-positive 4x4 matrix is a sum of a nonnegative matrix and a 
positive semidefinite matrix (see [5]), we know F = {(a, b) : A G £4,2}- 
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Figure 5: The pictures of the curve ip(a, 6) = and region F for co-positive matrices in 
Example 16.71 The left is for (i) , and the right for (ii) . 



(ii) Consider the symmetric matrices A parameterized as 



.4 



1 1 + a 1 + 6 1 + 6 1 + a 

1 + a 1 1 + a 1 + 6 1 + 6 

1+6 1+a 1 1+a 1+6 

1 + 6 1 + 6 1 + a 1 1 + a 

1 + a 1 + 6 1 + 6 1 + a 1 



When a = —2,6 = 0, it is the matrix associated to the Horn's copositive form (see Reznick 
[25]). The polynomial (p(a,b) defining equation (|6.3|> is 

a 9 • 6 10 • (a + 1) • (a + 2) 4 • (6 + 2) 4 • (2a 2 + 4a - 6) 5 (26 2 + 46 - a) 5 - 
(a 2 - 3a6 + 6 2 ) 7 • (6 2 + 26 - a)(2a + 26 + 5) • (a 2 + a6 + 2a + 6 2 + 26) 5 . 

The curves in the right picture of Figure [5] are defined by <p(a, 6) = 0. Let 



F 



(a, 6) G 



Aijxfx^ 

4<i,i<5 



is SOS in x 



It is an unbounded convex set. By the method used in Example 14.41 F is drawn in the 
shaded area of the right picture in Figure [5j Let G = {(a, 6) : A G £5,2}- Clearly, F C G and 
the boundary of G lies on <p(a,b) = 0. Then / a>6 (l, 1,0,0,0) > 0, f' a>b (l, 0, 1, 0, 0) > 0, and 
/ a ,&(l) 1,1,1,1) > imply that any pair (a, 6) G G satisfies 

2a + 26 + 5>0, a + 2 > 0, 6 + 2 > 0. 

Since (-0.5, -1.88), (-1.88, -0.5), (-1.3, -1.3) G (verified by GloptiPoly 3 [II]), from the 
right picture in Figure [5l we can observe that F is a maximal convex region that satisfies the 
above three linear constraints, excludes the previous 3 pairs and has the boundary lying on 
ip{a,b) = 0. So F = G. □ 
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7 Conclusions and discussions 



This paper studies the algebraic geometric properties of the boundary dPd{K). When 
K = W 1 , dPd(K) lies on an irreducible hypersurface defined by the discriminant of a single 
polynomial; when K is a real algebraic variety, the boundary dPd(K) lies on a hypersurface 
defined by the discriminant of several polynomials; when K is a general semialgebraic set, the 
boundary dPd(K) lies on a union of discriminantal hypersurfaces. General degree formulae 
for these hypersurfaces and discriminants are also proved. An interesting consequence of 
these results is that — log <p(f) can not be a barrier for the cone Pd{K) when ip(f) is required 
to be polynomial in /, but it would be a barrier if ip(f) is allowed to be semialgebraic. 

Given general multivariate polynomials /o, . . . ,f m , how to compute the discriminant of 
type A(/o, . . . , / m )? When m = 0, there are standard procedures for computing A(/o). 
However, to the best of the author's knowledge, this question is open for m > 0. In computing 
A(/) for a single polynomial /, it is typically non-practical to get a general formula for A(/), 
but if f(x) has a few terms and its coefficients have a few parameters, is there any practical 
method for evaluating A(/) efficiently? These questions are interesting future work. 

Acknowledgement The author would like very much to thank Bill Helton, Kristian Ranes- 
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References 

[1] D. Bertsekas. Nonlinear Programming, second edition. Athena Scientific, 1995. 

[2] F. Catanese, S. Ho§ten, A. Khetan and B. Sturmfels. The maximum likelihood degree. 
American Journal of Mathematics, 128 (2006) 671-697. 

[3] D. Cox, J. Little and D. O'Shea. Ideals, varieties, and algorithms. An introduction to 
computational algebraic geometry and commutative algebra. Third edition. Undergradu- 
ate Texts in Mathematics. Springer, New York, 1997. 

[4] D. Cox, J. Little and D. O'Shea. Using algebraic geometry. Graduate Texts in Mathe- 
matics, 185. Springer-Verlag, New York, 1998. 

[5] P. Diananda. On non-negative forms in real variables some or all of which are non- 
negative. Proc. Cambridge Philos. Soc, 58:17-25, 1962. 

[6] I. Gel'fand, M. Kapranov, and A. Zelevinsky. Discriminants, resultants, and multidi- 
mensional determinants. Mathematics: Theory & Applications, Birkhauser, 1994. 

[7] J. Gouveia, P.A. Parrilo, and R. Thomas. Theta Bodies for Polynomial Ideals. SIAM J. 
Optim., Vol 20, No.4, pp. 2097-2118, 2010. 

Preprint, 2009. larxiv: 0809. 34801 

[8] G.-M. Greuel, G. Pfister and H. Schoenemann. SINGULAR: A Computer Algebra System 
for Polynomial Computations. Department of Mathematics and Centre for Computer 
Algebra, University of Kaiserslautern. |http : //www . singular . uni-kl . de/ index . html 



29 



[9] J. Harris. Algebraic Geometry, A First Course. Springer Verlag, 1992. 

[10] D. Henrion. Semidefinite representation of convex hulls of rational varieties. LAAS-CNRS 
Research Report No. 09001, January 2009. l arXiv: 0901 . 18211 

[11] D. Henrion, J. Lasserre and J. Loefberg. GloptiPoly 

3: moments, optimization and semidefinite programming. 

http : //homepages . laas . f r/henrion/ sof tware/gloptipoly3/ 

[12] J. B. Lasserre. Global optimization with polynomials and the problem of moments. SIAM 
J. Optim., 11(3): 796-817, 2001. 

[13] J. Lofberg. YALMIP: a toolbox for modeling and optimization in Matlab. Proc. IEEE 
CACSD Symposium, Taiwan, 2004. [www . control .isy.liu.se/~j o hanl 

[14] E. Looijenga. Isolated singular points on complete intersections. London Mathematical 
Society Lecture Note Series, 77. Cambridge University Press, Cambridge, 1984. 

[15] J. Nie, J. Demmel and B. Sturmfels. Minimizing polynomials via sum of squares over 
the gradient ideal. Math. Prog., Series A, Vol. 106, No. 3, pp. 587-606, 2006. 

[16] J. Nie and K. Ranestad. Algebraic degree of polynomial optimization. SIAM J. Optim., 
20 (2009), no. 1, 485-502. 

[17] J. Nie, P. Parrilo and B. Sturmfels. Semidefinite Representation of the k-Ellipse. IMA 
Volume 146: Algorithms in Algebraic Geometry (Eds. A. Dickenstein, F.-O. Schreyer, 
and A. Sommese), pp. 117-132, Springer, New York, 2008. 

[18] J. Nie and B. Sturmfels. Matrix cubes parametrized by eigenvalues. SIAM Journal on 
Matrix Analysis and Applications, Vol. 31, No. 2, pp. 755-766, 2009. 

[19] J. Nie, K. Ranestad and B. Sturmfels. The algebraic degree of semidefinite programming. 
Mathematical Programming, Ser. A, Vol. 122, no. 2, pp. 379-405, 2010. 

[20] P. Parrilo. Semidefinite programming relaxations for semialgebraic problems. Math. 
Prog., Ser. B, Vol. 96, No.2, pp. 293-320, 2003. 

[21] P. A. Parrilo and B. Sturmfels. Minimizing polynomial functions. In S. Basu and L. 
Gonzalez- Vega, editors, Algorithmic and Quantitative Aspects of Real Algebraic Geom- 
etry in Mathematics and Computer Science, volume 60 of DIM ACS Series in Discrete 
Mathematics and Computer Science, pages 83-99. AMS, 2003. 

[22] M. Putinar. Positive polynomials on compact semi-algebraic sets, Ind. Univ. Math. J. 
42 (1993), 969-984. 

[23] K. Ranestad and H.C. Graf von Bothmer. A general formula for the algebraic degree in 
semidefinite programming. Bulletin of LMS, 41 (2009), no. 2, 193-197. 

[24] K. Ranestad and B. Sturmfels. On the convex hull of a space curve. Submitted to Doc- 
umenta Mathematica, arXiv : 0912 . 2986 



30 



[25] B. Reznick. Some concrete aspects of Hilberts 17th problem. Contemp. Math., Vol. 253, 
pp. 251-272. American Mathematical Society, 2000. 

[26] R. Sanyal. F. Sottile, and B. Sturmfels. Orbitopes. Submitted to Mathematika. 

[27] K. Schmiidgen. The K-moment problem for compact semialgebraic sets. Math. Ann. 289 
(1991), 203206. 

[28] B. Sturmfels. Solving systems of polynomial equations. CBMS Regional Conference Series 
in Mathematics, 97. American Mathematical Society, Providence, RI, 2002. 

[29] B. Sturmfels and C. Uhler. Multivariate Gaussians, semidefinite matrix completion, and 
convex algebraic geometry. To appear in Annals of the Institute of Statistical Mathemat- 
ics, featured section on Algebraic Methods in Computational Statistics. 



31 



